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Abstract. The entanglement dynamics of arrays of qubits is analysed in the presence 
of some general sources of noise and disorder. In particular, we consider linear 
chains of Josephson qubits in experimentally realistic conditions. Electromagnetic and 
other (spin or boson) fluctuations due to the background circuitry and surrounding 
substrate, finite temperature in the external environment, and disorder in the initial 
preparation and the control parameters are embedded into our model. We show that 
the amount of disorder that is typically present in current experiments does not affect 
the entanglement dynamics significantly, while the presence of noise can have a drastic 
influence on the generation and propagation of entanglement. We examine under which 
circumstances the system exhibits steady state entanglement for both short {N < 10) 
and long (N > 30) chains and show that, remarkably, there are parameter regimes 
where the steady state entanglement is strictly non-monotonic as a function of the 
noise strength. We also present optimized schemes for entanglement verification and 
quantification based on simple correlation measurements that are experimentally more 
economic than state tomography. 
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1. Introduction 

A fundamental property of the superconducting state is that it exhibits quantum 
coherence at the macroscopic scale, a feature that has been used to probe the validity 
range of quantum mechanics beyond the microscopic realm [U [2]. The development 
of quantum information science and the experimental progress in the manufacturing 
and control of superconducting-based quantum circuits has allowed for novel proposals 
aimed at implementing quantum information processing using Josephson qubits [3] . This 
generic denomination refers to qubit realizations that involve the charge [1] or the flux [5] 
degree of freedom in superconducting devices (also see References [HI [TJ El E]). The 
coherent coupling of two charge qubits and the implementation of conditional gate 
operations [10] , as well as the coupling of two flux qubits [11] , have been demonstrated 
experimentally, and there is currently an increasing activity in the field. Interesting 
applications include proposals to interface such devices with optical elements in order to 
create hybrid technologies [12]; or to use them for quantum communication [T3 | fH l [28]. 
It needs to be realised however that the technological barriers for full scale quantum 
computation are formidable. Thus there is a need for intermediate experiments that are 
interesting and non-trivial yet less demanding than implementing quantum computation. 
The exploration of many-body dynamics provides such a platform. Crucially, the 
fabrication of arrays that involve ~ 50 Josephson qubits has already been achieved 
in the laboratory [15], and one of our aims in the present work is to make realistic 
predictions about their dynamical entanglement properties. In order to do so, we will 
take into account the influence of certain general forms of noise and disorder on the 
system. 

A well- understood source of noise in any Josephson device is due to electromagnetic 
fluctuations in its background circuitry [16]. In the case of a single superconducting 
qubit, generic spin or boson fluctuations with a variety of spectral properties can be 
treated within the framework of the spin-boson model [1^. Note, however, that the 
precise sources of l//-type noise have yet to be identified and that the spin-boson 
formalism ceases to be valid in the limit of strong coupling to environmental fluctuators 
[T8] . Moreover, the influence of noise on A^ > 2 coupled Josephson qubits remains 
largely unexplored [19]. In relation to quantum information processing, it is important 
to characterise the necessary conditions for preserving coherence in a noisy environment 
before further steps can be taken in the direction of designing error correction schemes 
and (subsequently) fault tolerant superconducting architectures. 

In this paper we formulate an initial model for Josephson-qubit chains in realistic 
environments, taking into account the most common sources of noise. First we analyse 
the quantum dynamics in ideal conditions and then discuss the modifications one should 
expect when (i) disorder is taken into account and (ii) the system couples linearly to an 
environment that is modelled as a bath of harmonic oscillators, as detailed in section 4. 
To corroborate our view that the findings for shorter chains are generic, we also perform 
simulations for longer chains with A^ ~ 50 qubits. The simulations are performed 
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using a time-dependent Density Matrix Renormalization Group (DMRG) technique |2Uj . 
employing a code previously developed and tested in [21] . Given that our interest focuses 
on the study of quantum coherence, the system dynamics is characterized in terms of 
entanglement creation as well as entanglement propagation along the chain. There is 
currently no unique way to quantify entanglement in mixed states (see Reference [22] 
for a recent review) . We choose the logarithmic negativity [231 122] largely for its ease in 
computation and the availability of an operational interpretation [21]. It is defined as 

^N(Pi,j) = iog2 IIpJ^II , (1) 

where 1 1 . 1 1 denotes the trace norm of a matrix and pj'j is the partial transpose of the 
reduced density matrix pij for two subsystems i,j. Other measures would give broadly 
equivalent results [22] . 

Another fundamental problem concerns the development of techniques that allow 
for the detection and quantification of the entanglement that is present in a network 
of qubits. Exciting experimental progress in this direction for the case of Josephson 
qubits has been reported very recently [25] , whereby the entanglement was demonstrated 
via full state tomography. As the latter is a costly and time consuming experimental 
technique, strategies aimed at establishing a lower bound on entanglement by means of 
determining spin-spin correlations have been developed [26]. We test the performance of 
these concepts in the present case and find that, using some optimisation, they provide 
very accurate estimates for the amount of entanglement present in the system. 

2. Entanglement Dynamics under Ideal Conditions 

We consider an open chain of qubits with nearest-neighbour interactions. The 
Hamiltonian of the system is 

= -2 E (e.^f + A.^n - 2 E K^a! ® a^^, (2) 

1=1 i=l 

where cxf'^'^ denote Pauli matrices for qubit i, and Ki is the strength of the coupling 
between nearest neighbours i,i + 1 (we set /i = /c^ = 1 throughout). The control 
parameters are the energy bias and the tunnelling splitting Aj. We consider, as an 
example, charge qubits [3], in which case we have = 4:Ec{l — 2Ng) and Aj = Ej, where 
Ec is the charging energy, Ej is the Josephson energy, and 2eNg = CgVg is the gate 
charge, which is controlled by the gate capacitance Cg and voltage Vg. Charge qubits are 
operated in the regime where Ec ^ Ej\ therefore the energy scale is set by Eq-, which 
was of the order of 1 K in the experiment of Reference [4], and we let Ej/Ec = 0.1. 
We consider purely capacitive coupling between the charge qubits [9], and hence the 
(Ti ® o^^_^l interaction dominates. We assume (this condition will be relaxed later on) 
that the effective charge number of each qubit is Ng = 1/2 (i.e., = 0) so that it is 
operated at the so-called degeneracy point [7] . As it will become clear later, this choice 
can be advantageous when trying to minimise the impact of noise. 
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A feasible way to achieve generation of entanglement in coupled many-body systems 
is non-adiabatic switching of interactions as demonstrated in harmonic chains [271 128] . 
This approach has the advantage of only requiring moderate control over the parameters 
of individual subsystems. In our study here we will quantify the amount of entanglement 
that can be obtained in this way for the model Equation ([2]). In particular, we will 
assume that the interqubit couplings Ki are initially zero and are then non-adiabatically 
switched to a finite value. If one were indeed able to switch off the interqubit couplings 
completely, then at absolute zero temperature each qubit would be prepared in its 
ground state |+) (when operated at its optimal point, where e, = 0). The ground state 
of the Hamiltonian of Equation ([2]) for Ki = is the uncorrelated state 

l*(o)) = l+r^ i+) = ^(iT) + u)) (3) 

where | |), | |) denote the eigenstates of corresponding to zero or one extra Cooper 
pair in the superconducting box. We will thus study the generation of entanglement by 
evolving the initial state of Equation ([3]) according to the Hamiltonian ([2]) with Ki ^ 0. 
We will also study the propagation of entanglement [13l [HI [211 123 [28] by assuming that 
our initial state is 

\m) = i/5)i2 ® 1+)^''"', \p)i2 = ^(1 Ti) + 1 it))- (4) 

In this case, the interactions Ki are initially zero but a Bell state \(3)i2 has been created 
for the first two qubits. Again, the interactions are instantaneously switched on to a 
finite value and the evolution according to the Hamiltonian ([2]) with i^j 7^ is studied. 
The Bell state \P)i2J shared between the first two qubits in the chain, is maximally 
entangled. We note that there is also entanglement generation while the quantum 
correlations of \P)i2 propagate along the chain. We discuss below the effect of deviations 
in these initial conditions, due to non-vanishing initial interactions or static disorder. 

We begin by calculating the time-evolution of the logarithmic negativity of Equation 
([T]) for qubit pairs in ideal conditions. In Figure [1] we show the result for a chain of = 8 
qubits with the initial state |^I/(0)) of Equation ([3]) and parameters = 0,Ki = Aj/4. 
Due to the geometry of the setup, symmetric pairs of qubits, such as (1,2) and 
(7, 8) possess the same amount of entanglement. It is possible to create long-range 
entanglement, even between the first and last qubit in the chain (at t ~ 200 E^^ 
which corresponds to about 1.5 us). By comparing and contrasting the four panels in 
Figure 1 we can see that there is a characteristic 'collapse- and- revival' behaviour: when 
the entanglement of nearest or next-nearest neighbours is constant or vanishing, the 
entanglement of distant qubits becomes maximal (e.g., at t ~ 200 Eq^). Due to the 
finite length of the chain we find that neither the time when a pair of qubits becomes 
entangled for the first time, nor the time when the first entanglement maximum occurs, 
are proportional to the distance between the qubits. The characteristic speed at which 
the distance of pairwise entanglement is expected to grow is thus masked by finite size 
effects, in our study. Entanglement propagation for this chain in ideal conditions is 
considered later ( cf. Figure [5]) . 
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Deviations in initial conditions: In practice, it is not quite possible to switch off 
the interqubit couphngs completely. To take this into account we have considered the 
case when there is initially some small coupling between the qubits, -ft^mi, and the initial 
state of the system is the ground state of ?is(-f^ini)- Then the ground state evolves 
according to 7is(-f^fin), where Kfi^ = A/4. In this case we obtain very similar results 
with those presented in Figure [1] for the ideal case (clearly, for Ki^^ — > we recover the 
results of the ideal case). In particular, for A'ini < A/100 the relative maxima deviate 
by less than 5%, and there is initially very little entanglement in the system (e.g., the 
logarithmic negativity for the first two qubits in the chain at t = is less than 0.004). 
By contrast, for i^ini ~ A/10 the relative maxima can deviate by up to 30% and the 
initial entanglement in the chain is much more evident (e.g., E^{pi 2) ~ 0.1 at t = 
for the same parameters). We will revisit this case shortly, after we have introduced 
disorder and noise into the system (c/. Figure HI). 

Another interesting question relating to the initial preparation concerns the state of 
the chain at thermal equilibrium. In particular, we would like to know if we would obtain 
similar results when the state of the system at t = is the thermal equilibrium state, 
and also how close are the thermal and ground states of the system described by the 
Hamiltonian 7Ys(-^ini) of Equation ([2]). We have therefore assumed that the initial state 
of the chain is the thermal equilibrium state p(T) = exp{— /3T-Cs) / Z , where /3 = 1/T, 
Z = Tr[exp(— /37is)], and there is initial coupling, Ki^i, between the qubits. The coupling 
is then switched on to its final value i^gn at t = 0. In this case we have found that for 
low temperatures (T < 20 mK) the entanglement dynamics of the chain is very similar 
with that obtained by evolving the ground state (the relative deviations are less than 
10%) for the same value of i^ini- In order to compare the thermal equilibrium state 
p(T) with the ground state \G) of HsiKi^ii) we have calculated the fidelity {G\p(T)\G) 
for various values of the temperature (with fixed K[^i = A/4). We have found that for 
temperatures T < 15 mK the fidelity is between 0.99 and 1, and hence the two states 
are very close indeed for these temperatures. Between 15 mK and 25 mK the thermal 
equilibrium state and the ground state begin to differ (their fidelity slowly drops to 0.9 
as the temperature is increased). 

3. Influence of Disorder 

In any experimental situation the initial preparation will also suffer from errors in the 
control parameters Octri = Ci, Aj, i^j. As a result, the homogeneity of the chain will 
be broken. We can simulate the effect by letting the parameters take random, but 
static, values in the interval [(1 — d)actrh (1 + d)actr\], where d quantifies the disorder. 
An example is shown in Figure [2](a), where we plot E^{t) for the pair (1,2) in the 
ideal (solid line) and imperfect scenario (broken line), where the disorder in Aj and 
Ki is d = 0.05. Averaging over 10^ runs, we have found that disorder with d = 0.01, 
d = 0.05, and d = 0.10 causes relative fluctuations of the maximal entanglement equal 
to 0.011, 0.027 and 0.054, respectively. Therefore for disorder which is less than 10% 
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Figure 1. Pairwise entanglement creation for a homogeneous chain of iV = 8 qubits in 
ideal conditions; the system is described by the Hamiltonian of Equation ^ and the 
initial state is |5'(0)) of Equation ([3]). Note that symmetric qubit pairs, such as (1,2) 
and (7,8), generate the same entanglement and hence are represented by the same line. 
A 'coUapse-and-revival' pattern emerges in the entanglement oscillations, as seen from 
a comparison of the different panels. 



(the upper bound in state-of-the-art experiments [29]) the entanglement in the system 
changes marginally, on average. This is indeed true for the noisy scenario also, as shown 
in the following section. It is noted that disorder has recently been studied in related, 
but different, contexts in [271 [30] . 

4. Noise Model For a Variety of Sources 

We consider a spin-boson Hamiltonian of the form 

N 

'HsB = 'Hs + 'Hb + ^ crtXi (5) 

i=l 
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Figure 2. Entanglement creation between qubits (1, 2) in the presence of (a) disorder 
and (b) noise, which is characterised by a decay rate F at absolute zero temperature. 
Subplot (c) shows E]\f{t) for qubits (1,8) for various values of the decay rate F, 
temperature, and disorder. In (a) the initial state is |^(0)) of Equation ([3]) and it 
evolves under the Hamiltonian ([2]); in (b) and (c) the initial state is IV'(O)) of Equation 
(|5]) and it evolves under the master equation ©. 



where the first term corresponds to the free system Hamiltonian of Equation ([2]), 
the second term is the Hamiltonian for all independent baths i = {1, 2, . . . , A^}, 
Hb = J2iLiJ2k^k^(^k^^(^k\ where the k-th mode of bath i has boson creation and 
annihilation operators a^*^^ and a^^\ respectively, and the third term is the interaction 



between a qubit and its bath, whose 'force' operator is Xi = Y^iQi aj + ct/ [T7] . 

5- 



Clearly, it is assumed that each qubit is affected by its own bath, i.e., 0-^ ■'^] = ^ij 



a reasonable requirement for charge qubits biased by independent voltage gates. 

In the coherent regime, where uJi = (ef + A^)^/^ is much larger than the thermal 
energy, the preferred basis is given by the eigenstates of the single-qubit Hamiltonian, 
i.e., |0) = cos(^i/2)| t) + sm{e,/2)\ [) and |1) = -sin(^i/2)| t) + cos(^,/2)| i), where 
the mixing angle obeys tan 6*4 = Aj/ej. In this basis, TisB becomes 

N 

Kb = K + + E(sin + cos e^a^)X, (6) 

i=l 

where 

N ^ N-1 

= -9 E ^i^' - 9 E K^ic^Crt + Si(Tf )(Ci+iCr^Vl + S.+ldf+i) (7) 
^ i=l ^ i=l 

is the system Hamiltonian (the Pauli matrices are now written in the {|0), |1)} basis) 
with Cj = cosftj. Si = sinftj. In this basis the states |^I/(0)) and |$(0)) of Equations ([3]) 
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and (jl]), respectively, become 

\m) = \Of\ |0(O)) = |/3')i2®|O)^^-' (8) 

where \(3')^^ = + |10)). 

When the bath's degrees of freedom are traced out, and within the Born-Markov 
approximation, the time evolution of the chain is governed by a master equation of the 
Lindblad form 

p = -i[H',,p]+Cp (9) 

where is given by the same expression of Equation (7) provided that the weak 
coupling limit, where Ki < Ui, holds. The damping terms are given by the usual 
expressions, 

TV 
i=l 

+gi{2atpat - 2p)] (10) 

where af = 2~^{a^ ± zaj) and the parameters are defined as 

Gi = sin^ 9i{l + nT)r, Gj = sin^ OiUTT, gi = cos^ 9iT (11) 

with nx denoting the average number of bosons in the environment. We assume that the 
environments of all qubits are identical. We have not specified the environment's spectral 
properties and hence the decay rate T is given as a phenomenological parameter whose 
exact value can be adjusted to match that obtained for the actual spectral density of the 
bath. Within this framework, a broad class of dissipative effects can be accounted for, 
ranging from electromagnetic fluctuations in the surrounding circuitry to, for instance, 
the coupling to a phonon bath [31]. The system can therefore be viewed as a non 
critical dissipative Ising chain (The critical behaviour of dissipative Ising chains has 
been recently addressed in [32].) 

At the degeneracy point ej = and cos6'j = 0. As a result, each qubit is susceptible 
to relaxation only (the 'optimal' point introduced in Reference [7]). If, however, the 
energy bias is not exactly zero then the longitudinal contribution o^^Xi leads to pure 
dephasing at a rate g^. In any experimental realization, the presence of disorder limits 
the accuracy with which qubits can be operated at their optimal points, especially when 
it comes to the operation of long chains. In what follows we take this into account and 
study the modifications due to disorder. In current experiments [29] the value of disorder 
is typically 5 — 10% at temperatures 20 — 40 mK. The decoherence time td = T"^ for two 
coupled charged qubits [10] was reported to be around 2.5 ns (for single qubits td can 
be higher). In our simulations below we usually assume a worst-case scenario and let 
d = 5% at T ^ 41 mK with decay rate T = lO^^E'c (which corresponds to td ~ 1 '^'S). 

5. Dynamics of Short Chains (A^ ~ 10 qubits) 

In this section we discuss the results on the entanglement dynamics of open chains 
with N = 8 qubits (Figures [TH). Figures [1] and [2]|^a) have been analysed previously. 
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Figure 3. Entanglement creation between two qubits for the case of relaxation only 
(top) and relaxation, finite temperature, and disorder (bottom). The initial state is 
|^/'(0)) of Equation ^ and the system evolves under the master equation 



Figure[2](b) shows the creation of entanglement in the pair (1, 2) for various values of the 
relaxation rate F. Figure [2](c) shows the creation of entanglement in the pair (1,8) for 
various values of F and other parameters. It is seen that for noise strength F = 5 x 10~^ 
(i.e., td ~ 5 ns), average number of photons = 0.01 (i.e., T ^ 22 mK) and disorder 
d = 5% one may still obtain substantial entanglement between the first and last qubit 
in the chain (in particular, the ratio of the values of the first maxima corresponding to 
the imperfect / ideal cases is approximately 2/5). 

In Figure [3] we plot E]y(t) for different pairs of qubits in the case of (top) relaxation 
with F = 10-2 at zero temperature (T = 0) and (bottom) relaxation with T = 10 ^, 
finite temperature T ^ 41 mK (tt-t = 0.1) or T ^ 33 mK [ut = 0.05), and disorder 
5% {d = 0.05). As expected, the entanglement beyond nearest neighbours is drastically 
reduced in the presence of larger values of the noise (the correlations between the first 
and last qubit in the chain vanish altogether for this high value of F). 

In Figure H] we study the creation of entanglement between the first two qubits in 
the chain (in subplots (a) and (b)) when the initial state is the ground state of the 
Hamiltonian 7ig(-R'ini) of Equation ([7j), for various values of the initial coupling strength 
Ki^i. At t = the coupling is instantaneously switched on to its final value Ka^ = A/4. 
Subplot (a) shows the case with noise, at absolute zero temperature. Subplot (b) takes 
into account the temperature in the environment (T ^ 41 mK) . In subplot (c) we study 
the case whereby at t = the state of the system is in thermal equilibrium with its 
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Figure 4. Entanglement creation between qubits (1,2) in noisy conditions when 
there is some initial homogeneous coupling Ki^i that is instantaneously switched on 
to -ftTfin = A/4 at t = 0. In subplots (a) and (b) the initial state is the ground state 
('GS') of Ti,'g{Kini) of Equation ([7]) and they correspond to temperatures T = and 
T « 41 mK (riT = 0.1). In subplot (c) the initial state is the thermal equilibrium state 
('TES') of the system at a given temperature T, for various values of the parameters. 



environment at temperature Tq. Therefore in this case we let p(To) = exp(— /57ig)/Z 
at t = 0. The system Hamiltonian Tig, given by Equation ([7]), depends on the initial 
interqubit coupling Kini- The evolution proceeds according to the master equation 
(Q with the coupling Kf[^ and an average number of photons that corresponds 
to the temperature Tq. It is seen that at operating temperatures of around 40 mK 
the entanglement vanishes. It is however possible to observe entanglement when the 
temperature gets smaller (e.g., for Tq 33 mK). 

The results in Figure H] seem to indicate that the increase in the noise strength 
r and the external temperature yield the unavoidable degrading of entanglement 
generation. The amplitude of the entanglement oscillations decreases and the system 
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Figure 5. Steady state entanglement between qubits (1,2) in a chain of iV = 4 when 
riT = 0.1 and for different values of the ratio K/A as a function of the noise strength 
r. Parameter regimes can be identified where entanglement generation is enhanced by 
amplifying the environmental noise. 

becomes separable in the steady state for F and/or nx sufficiently large. However, this 
behaviour is not universal and we need to differentiate two time scales in the system. 
The initial transient is always such that the amplitude of entanglement oscillations is 
reduced as the noise increases and the amplitude of the first entanglement maximum 
is a monotonically decreasing function of both F and T. However, for a fixed riT, the 
steady-state entanglement can display a non-monotonic behaviour as a function of F. 
This phenomenon is of the same type of the noise-assisted effects that have been studied 
in Reference [33] for weakly driven spin chains and is illustrated in Figure [5] for a system 
of = 4 qubits. We see that at the selected temperature where rix = 0.1, there are 
parameter regimes for which the steady-state entanglement is initially zero for low values 
of the noise strength and resurfaces when F is increased over a certain threshold. This 
result indicates that if the aim is to generate entanglement in the steady state, it may 
be advantageous to amplify the environmental noise so as to maximise entanglement 
production along the chain. Persistence of this effect in longer chains (A^ ~ 40) has 
been corroborated numerically. 

Propagation of entanglement is analysed in Figure [6] for (a) ideal and (b) non- 
ideal conditions. In the ideal case, entanglement propagates from the ffist two qubits 
to the last two, but not perfectly. When we take into account noise and disorder 
the entanglement transfer is not possible and the last two qubits quickly reach their 
steady state, which is slightly entangled at absolute zero temperature, but separable at 
T ~ 20 mK for the selected parameter regime. 
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Figure 6. Entanglement propagation at zero temperature for the ideal case (top) 
and for relaxation and disorder (bottom). In the top plot, the initial state is |5'(0)) 
of Equation (|4]) and the system is described by the Hamiltonian Hs of Equation ([2]). 
In the bottom plot the initial state is |(?I>(0)) of Equation ([8]) and it evolves under the 
master equation ([9]). 



6. Dynamics of Long Chains 

To confirm the vahdity of our findings for longer chains, we have performed time- 
dependent DMRG simulations [20] . 

For an ideal chain without noise and disorder, we have considered entanglement 
generation in the model ([2]) with = 20 qubits. Here the matrix dimension was chosen 
dim = 20 and a 4th order Suzuki- Trotter decomposition was employed. The results 
are in good agreement with the findings for shorter chains in figure [H In particular, we 
observe the same 'collapse-and-revivaF pattern of the entanglement oscillations, and the 
long-range entanglement peaks at those regions where the short-range entanglement is 
close to its steady-state value or vanishes altogether. 

For the cases which include noise and disorder, a matrix product representation 
for mixed states with matrix dimension dim = 60 and a 4th order Suzuki- Trotter 
decomposition were used for a chain of = 40 qubits [20] . A sketch of the method is 



given in Appendix A 



Figure [7] shows the creation of entanglement in the presence of noise, at zero 
temperature for both a homogeneous and a disordered chain (in which case disorder 
occurs in as well as Ai,Ki). Figure [HI shows entanglement creation in a noisy 
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N = 40; A' - 0.025;& 0.01;nr = 0; dis = 0% N ^ 40; A' ^ 0.()25;I^ O.Olmr ^ 0; dis ^ 5% 
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Figure 7. Entanglement creation, i.e. evolution of the initial state |V'(0)), given 
by Equation ([9]) in a chain of TV = 40 qubits at zero temperature, in the presence of 
noise (F — 0.01). The two plots on the left show the homogeneous case, while the two 
plots on the right show a case with 5% disorder in ei, Ai, Ki. Qubits at the boundaries 
are slightly stronger entangled than in the centre of the chain. The entanglement 
between qubits that are further apart than shown here is zero. 

homogeneous chain for various values of temperature. For all quantities we find good 
agreement with the results obtained for N = 8, where the relative deviations between 
N = 40 and N = 8 are less than 5%. It is also noted that the entanglement between 
two blocks of two qubits each was found to be about 17% higher than the entanglement 
between individual qubits of the same separation. 

7. Witnessing Quantum Correlations: Experimental Verification of 
Entanglement 

In experiments it will be crucial to verify the existence of entanglement via 
measurements, which ideally should also permit a quantification of the detected 
entanglement. This could be done by full state tomography, which is a very costly 
experimental procedure though. Being able to establish a lower bound on entanglement 
from the measurement of a few observables will thus be a significant advantage. Recently, 
a theoretical framework for the exploration of these questions has been developed for 
general observables [26] and witness operators |34j . 

The basic approach is to identify the least entangled quantum state that is 
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Figure 8. Entanglement creation, i.e. evolution of the initial state \ip{0)), ([8]), given by 
equation ([9]) in a chain of = 40 qubits at various temperatures (nx = 0.05; = 0.1 
and riT = 0.2) in the presence of noise (F = 0.01). Only nearest neighbours become 
entangled in this case. 



compatible with the measurement data. The entanglement of that state then provides a 
quantitative value for the entanglement that can be guaranteed given the measurement 
data. In [26], in particular, spin-spin correlations have been used to determine such a 
lower bound analytically. We now employ this concept for our system and consider the 
two quantities, 

Ci(p„-) ^ max [0,log2 + |q,^|)] (12) 



= max [0,log2 (l + iqjl + l^f | + - ij (13) 

where = Tr[afa'^p] {a,b = x,y,z). Both quantities form a lower bound to the 
logarithmic negativity, i.e. E^{pi j) > Ci(pjj) and E^{pi j) > C2{pij). Ci{pij) can be 
employed if only Cfj and C?j are accessible in measurements. However, if Cfj can be 
measured too, then C2{pij) yields a tighter bound. 

Figure [9] shows that both lower bounds provide good approximations for the 
logarithmic negativity of two neighbouring qubits. If the qubits are next-nearest 
neighbours, C2{pij) still provides a good estimate, while Ci(pjj) eventually fails to 
approximate the entanglement well. 

The reason why Ci{pij) and C2{pij) sometimes do not approximate the 
entanglement very well lies in the choices of the axes along which correlations are 
measured. Instead of Cfj, Cfj and C^j one could consider correlations along a rotated 
set of axes, C^^, C^j and C^j, where = J2a=x,y,z ^aaO'T and Raa is an orthogonal 
matrix representing the rotation. Choosing to measure correlations along x, y and z 
may hence underestimate the entanglement severely. The best approximation of the 
entanglement is obtained by maximising Cf", C^^j and Cf^ over all possible choices for 
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the axes a, b and c. 

This optimal choice can be obtained in the following way: If the state p is symmetric 
with respect to subsystems i and j in the sense that C^^J = Cfj, Cfj = Cf^ and 
Cfj = C^j, then the matrix 



X 



I /^xx r^xy /^xz \ 

I \ 

^yx ^yy ^yz 

i,j i,j i,j 



(14) 

is real and symmetric and hence has real eigenvalues and is diagonalised by a rotation. 
Let us denote the eigenvalues of X by Ai, A2 and A3, then the quantity 

C^(pi,,-) =max[0,log2(l + |Ai| + |A2| + |A3|)-l] (15) 

provides the best approximation of E^{pij) of the form (fT2l) . as Ai, A2 and A3 are the 
spin-spin correlations along the optimal choice of axes 

As an example, in figure [9], the entanglement between qubits 10 and 11 is 
^n(Pio,ii) = 0.2096 at t = 10E^\ While C2(pio,ii) = 0.1583 at this point, we obtained 
C2(pio,ii) = 0.2096 for the optimal choice of axes a, b and c. The optimal choice of axes 
depends on time. Yet one fixed set of axes approximated the entanglement very well 
over a range of At = 5Eq^ in our example. 

8. Conclusions 

We have studied the dynamics of entanglement in qubit chains influenced by noise 
and static disorder. This study provides useful analysis for interesting experiments 

I According to theorem VIII. 3. 9 of [35j, for a square matrix X, J^i l^ul ^ Si I'^il; where the Ai are the 
eigenvalues of X. As |Ai| is a unitarily invariant matrix norm, it does not depend on the choice of 
basis while J2i l^ul does. Thus the largest value that can be achieved for J^i is given by J2i 
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with quantum devices that in the long term may be suitable for the implementation of 
quantum computing. We have considered an experimentally interesting implementation 
using Josephson charge qubits with capacitive interactions between nearest neighbours. 
We have found that static disorder less than 10% (i.e., the current experimental upper 
bound) does not affect the entanglement dynamics substantially. By contrast, the 
influence of environmental noise, modelled here as a set of independent harmonic 
oscillator baths of arbitrary spectral density, is much more pronounced: it reduces long- 
range correlations and decreases the magnitude of the achievable bipartite entanglement. 
For typical operating temperatures, the influence of noise on the chain dynamics at 
short times and in the steady-state can be crucially distinct; while the entanglement 
amplitudes in the initial transient decrease monotonically with the noise strength, the 
steady-state response is non-monotonic. Therefore, we have identified parameter regimes 
in which the bipartite entanglement increases as a result of amplifying the noise. We 
have found agreement between the behaviour of entanglement in short (A^ ~ 10) and 
long (A^ ~ 50) chains. 

The present results are encouraging from the experimental point of view as 
they suggest that (a) both short- and long-range entanglement can be generated and 
propagated under suitable laboratory conditions, and (b) lower bounds can be placed 
on the entanglement from the measurement of a few general observables. 
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Appendix A. Matrix Product State Simulations for Mixed States 

Here we outline the concept proposed in [20] and its adaption to our application. For 
the Matrix Product State simulation of mixed state dynamics, the density matrix for A^ 
qubits is expanded in a basis of matrices formed by direct products of the elementary 
matrices 




(A.l) 




(A.2) 



Hence the matrices forming the basis for A^ qubits are of the form 



® (g) . . . (g) . 



(A.3) 



A*" sites 
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The expansion of the density matrix p is now written in terms of products of matrices 
in the following way: 

P= E ■ Ai ■ rf ^1 ■ A2 ■ . . . ■ A^_i ■ r[^"l e,, ® e,, ® . . . ® e,^ (A.4) 

si,S2,.--sjv=l 

where '■' denotes matrix multiplication. Here, each r[*^' (si = 1,2,3,4) is a row vector 
of length D, each fJ^'^' is a column vector of length D, each F^*^' (j 7^ 1, A^) is a, D x D 
matrix and each Aj is a diagonal D x D matrix. The structure of the matrices F and 
A is the same as in the Matrix Product representation of pure states and the TEBD- 
algorithm [20] can be employed for the simulation of the dynamics. In contrast to 
pure states, the matrix elements of the Aj for mixed states can however no longer be 
interpreted as the Schmidt coefficients of the respective decomposition. 
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